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We demonstrate the appearance of unexpected reflection resonances in corrugated dielectric 
waveguides. These are due to the curvature of the boundary. The effect is as strong as the or¬ 
dinary Bragg resonances, and reduces the transmission through our waveguide by 20%. It is thus 
of high relevance for the design of optimized waveguiding structures. We validate our analytical 
predictions based on coupled mode theory by a comparison to numerical simulations. 
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I. INTRODUCTION 

Light scattering is the key mechanism to tailor the 
properties of passive waveguiding structures. Ordered 
and disordered boundaries determine the transmission, 
the reflection, and the radiation losses. In highly in¬ 
tegrated structures these processes can be controlled 
through designed lithography, to modify light propaga¬ 
tion at will. For this reason, a thorough understanding 
of scattering processes is mandatory. One powerful tool 
for understanding the light propagation and scattering 
is the concept of Bragg scattering. Given only the peri¬ 
odicity of a waveguiding structure, it predicts reflection 
resonances in a straightforward manner. Therefore, light 
propagation through dielectric waveguides, has been in¬ 
tensively investigated concerning Bragg scattering m- 

In this work we show that care has to be taken when 
applying Bragg scattering analysis in a too simplified 
way. Surprisingly, rather strong reflection resonances can 
be overseen. 

In a waveguide with periodic boundaries the m-th or¬ 
der Bragg reflection is in general expected for the wave¬ 
length A that (for perpendicular incident) fulfills 

2d = m\ (1) 

where d is the length of the periodicity, e.g., the lattice 
spacing. 

The appearance of multiple order m seems obvious, be¬ 
cause reflection takes place, whenever the backscattered 
wave interferes constructively with the incident wave. 
This is always fulfilled for wavelength increments of 2d. 
Surprisingly this simple picture of multiple order is in¬ 
correct for some systems. There are periodic systems 
where only a single Bragg reflection (m = I) exists. For 
example, waveguides with infinite sinusoidal boundaries 
(similar to the finite waveguide sketched in Fig. [^. The 
Fourier series of the sinusoidal boundary consists of only 
a single (positive) coefficient. Therefore these systems 
are believed to show a single Bragg reflection only [7]. 
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Figure 1. Model of corrugated waveguide. Transmittance T, 
reflectivity R, and losses l — {R-\-T). The (mean) width of the 
waveguide is d = 450nm, the length L — lOd = 4.5pim. The 
refractive index of the inner (outer) material is n = 2 (n = 1). 
The wavelength of the boundary oscillation is A = 200nm, 
and its amplitude is A = 37.5nm {a = 37.5nm/%/2). 


Here we perform a more rigorous analysis and show 
that this is in fact not true. It turns out that corrugated 
waveguides show an additional Bragg reflection which is 
not expected from previous studies on individual peri¬ 
odic systems [ZllHl. Since the simple Bragg picture in 
Eq. 0 often serves as starting point for numerical de¬ 
sign of optical components Eli, it is vital to know which 
resonances can exist in principle. 

Fig. [2 shows a typical finite waveguide which we study 
in our analysis more explicitly as an example. We apply 
a technique from statistical boundary roughness analy¬ 
sis [3] to single dielectric waveguides. Within a coupled 
mode approach we are able to generalize findings for en¬ 
sembles of systems to single systems under drastically 
relaxed assumptions. The previous theoretical findings 
become a special limiting case in our framework. 

The evidence of an additional Bragg resonance pre¬ 
sented here, is important for many systems in a variety of 
communities, where corrugated waveguides are employed 
in very different applications, such as group-velocity con¬ 
trol [1], phase-matching in nonlinear materials [2], dis¬ 
tributed feedback laser [4], optical filtering [3], grating 
couplers [S], and hybrid atom-photonic systems m- 
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II. COUPLED MODE APPROACH 

Coupled mode theory is a powerful and common 
method j7]. The starting point is the wave equation for 
a dielectric waveguide with permeability /i and dielectric 
function e(a;, y), which is weakly perturbed by As{x, y, z), 
written as 

[\/'^ + uj'^n{e + Ae)] E{x,y,z) = 0 . ( 2 ) 

The field E of frequency ui = kc can be constructed 
from the Eigenfunctions Ei of the unperturbed waveg¬ 
uide. The contribution Ai of each Eigenmode Ei , changes 
along the waveguide, due to the dielectric perturbation, 
such that 

E = ( 3 ) 

i 

where /3m is the propagation constant of the m’th mode. 
This is, the component of the wave vector in propagation 
z-direction. Inserting Eq. (|^ into Eq. (H and taking the 
scalar product with E^ yields (see [7] for details) 

where the changes in Ai are assumed to be sufficiently 
“slow” to neglect (see [7]). The coupling coefficient 


where Cki is the Fourier coefficient in the Fourier series 
of Cm (as in Eq. ®)- 

A common simplification is to assume that the dielec¬ 
tric perturbation is periodic in z 0151 

Ae(a;,y,z) = 5]e„(x,j/)e-™T^ 

m 

Then the z-dependence of Ae{x,y,z) can be separated 
from the overlap integral in Eq. which is only over x 
and y. 

^ / dxdyEkemix,y)Ei ( 10 ) 

This means, in particular, that the integral becomes in¬ 
dependent of z. In fact we reduced the description of 
coupling effects to a stratified waveguide. This is, a 
waveguide which is composed of rectangular slices, each 
with some Ae. This is not surprising, because we con¬ 
structed the E-Field as contributions of Eigenfunctions, 
weighted by Ai in Eq. (§. We assumed that the {x, y, z)- 
dependencies of the E-Field can be separated into A;(z) 
and the Eigenmodes Ei{x,y). However, if we then try 
to tackle an arbitrary corrugated waveguide, we have to 
keep in mind that the only results we can expect are 
results for stratified waveguides. Common text book ap¬ 
proaches ignore this fact in lain]. 

III. COORDINATE TRANSFORMATION 


Cki j dxdyEkAeEi ( 5 ) 


describing the overlap between two modes k and 1. The 
coefficient u accounts for different possible choices of the 
normalization of E^. In the following, we will restrict 
ourselves to two modes. In this case Eq. Q is a set of 
two coupled equations (fc = 1, Z = 2 and fc -O- Z) 
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( 6 ) 

(7) 


A relevant figure is the reflectivity of a system. Solv¬ 
ing these coupled mode equations yields two solutions 
Ai,A 2 . The ratio of these two solutions at the begin¬ 
ning of the sample (z = 0), is the ratio of incoming to 
backscattered mode, i.e., the reflectivity 


A2(0) 

Ai(0) 


( 8 ) 


When solving Eq. §) and Eq. 0 and plugging the so¬ 
lution into Eq. Q, the maximum reflectivity is obtained 
as 


R 


max 


= tanh^ 



(9) 


Corrugated waveguides feature rich physical effects 
which are not present in stratified waveguides. There¬ 
fore, as a next step we will show how to overcome the 
shortcomings of common coupled mode theory, without 
rejecting the entire approach, which, indeed, captivates 
by its clarity and simplicity. 

As the above problems are of geometrical nature, i.e., 
restriction to stratified waveguides, it is reasonable to 
look for a geometrical solution. Here we use a straight¬ 
forward coordinate transformation to transform the cor¬ 
rugated boundaries to flat boundaries. 

This transformation has been used previously to derive 
the square gradient scattering mechanism in systems that 
feature boundaries with randomized roughness [9j. Even 
though this was the very first derivation of the square 
gradient mechanism, the general validity of this mecha¬ 
nism was so far doubtful for several reasons. At first it 
is derived for ensembles of systems that feature peculiar 
statistical properties. Therefore it is not a priori clear if 
the mechanism is only a statistical effect. Previous ex¬ 
periments [12] investigated systems that resembled these 
special statistical properties. Furthermore, both theory 
and experiment have so far being restricted to hollow 
waveguides with perfectly electrically conducting bound¬ 
aries. 

In the following, we will put the square gradient scat¬ 
tering mechanism on solid theoretical grounds for indi¬ 
vidual systems, with arbitrary boundaries. This will be 
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done for dielectric waveguides, but we will show that our 
general results are valid for perfect electric conductors as 
well. 

We assume a waveguide as shown in Fig.[2 The waveg¬ 
uide has a width d whose boundaries are given by a nor¬ 
malized boundary function q{z), such that the bound¬ 
aries are at a: = ±dl2 ± aq, where tr^ is the variance of 
the boundary. For a sinusoidal boundary a is connected 
to the amplitude A of the oscillation of the boundary by 
a = Ajspi. We chose the coordinate transformation as 

UJ (z) 

{x,y,z) —» ( ^ x,y,z) ^\th w{z) = d + 2aq{z) (11) 


It will flatten the boundaries of the dielectric waveguide 
at a; = ±d/2 and thus set Ae = 0 in Eq. ([^, yielding 


-b wVeo(i,y) E(x,y,z) = 0 


( 12 ) 


The transformed Laplacian consists of several new 
terms 


= V 


red f -b -b Vgg 


(13) 


which are calculated in Appendix A Here, the reduced 
Laplacian is defined as 


^red • — 


d3A dy"^ 


dP 


(14) 


The other terms will now be interpreted as the new di¬ 
electric perturbation Ae 


V? + 


=: w nAi 


(15) 


yielding the transformed wave equation 

Vred + (eo(i, y) + Af)] E{i, y,z) = 0 (16) 


In the next three sections the three terms in Eq. (13) 


will be studied in detail. It will be shown that yields 
the well-known Bragg reflection (therefore index b). It is 
analog to previous stratified approximations. Erequency 
analysis shows that can be safely neglected. Einally 
the last term, will turn out to represent the novel 
mechanism of square gradient Bragg reflection (index sg). 


IV. STRATIFIED APPROXIMATION YIELDS 
BRAGG SGATTERING 


To work out the physical influence of we ignore the 
two other terms in Eq. (15) and set 


A'= — 
tdV 


dP 


The coupling coefficient from Eq. becomes 


^O) _ - 
kl ^ 


dxdy—EkAeEi 

d 


(18) 


(19) 


where the prefactor is the Jacobian dxdy = dxdy | ^ |, 
which is strictly positive, so that |^| = The modes 

of the electric field are the undisturbed modes of 


the transformed system, i.e., Ae = 0 in Eq. (161. These 
modes are calculated in [Appendix B| 

Plugging Eq. (|I^ into Eq. (191, yields 


I didyE^i^E ( 20 ) 


puj^ 




Approximating 


dP 


a d + aq 




we have 




A) 

kl 


( 21 ) 


where 


Ab) _ g 4 

— 7 

a pujp 


- ^2 ^ 
dxdy Ek—El 


( 22 ) 


The coupling coefficient can now be readily calculated 
for given material parameters. Before doing so, we show, 
that the first term in Eq. (151 indeed corresponds to 


Bragg scattering. To this end, we analyze the period¬ 
icity of , which is clearly the same as the periodicity 
of the corrugated boundary q{z). Assuming that q{z) is 
a periodic function it can be expanded as 


qiz) = 




(23) 


We can thus write 
Ab) 




(24) 

(25) 


The first term 


V? 



(17) 


and set = qm I■ 

As in [7], we now analyze small changes in the ampli¬ 
tude A[z) in Eq. ([^ by integrating over a length s, which 
is long compared to A 


contains no derivative of q{z). In this term the curvature 
of the boundary has no influence. The physical effects 
expected from are those of the stratified waveguide. 


dz (26) 

l m ds 
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This integral will vanish unless the Bragg condition 

Pk- I3i = m— (27) 

is satisfied. In case of backscattering into the same mode 
(/3fc = —/3;)the Bragg condition takes the form 

Pk = (28) 

which is Eq. Q for lattice spacing d = A. 

We have thus shown, that the first term of the trans¬ 
formation yields the well known Bragg scattering. What 
about the different order ml For a sinusoidal bound¬ 
ary q{z) ~ sin( ^z) , there are only two possible values 
m = ±1 in Eq. ( |23[ ), which yield the same wavelength. 
Consequently, there is only one single reflection resonance 
predicted in the stratified approximation. 

After we have successfully recovered the results for the 
stratified approximation we will now turn to the next 
term. 


V. BEYOND BRAGG SGATTERING 


The second term V^, comprises derivatives of q{z). 


w \dz 


JLJL 

dx dz 


d^q^ d 


5 

dz dx 


Here the curvature of the boundary influences the scat¬ 
tering. This is the first indication that we no longer deal 
with a stratified system. However terms in will have 
the same periodicity as q{z) itself and come into play at 
frequencies given by the Bragg condition. This means 
that they act at the exact same frequencies as the terms 
used in the stratified approximation. We will neglect this 
term here, because in this study we are interested in re¬ 
flection resonance that take place at frequencies other 
than the Bragg resonances. 

In contrast, the third term on the left hand side of 


Eq. (15) contains terms with the square of the derivative 
of q{^. In general the square of a function can have 
a different periodicity than q{z). Therefore we expect 
this term to play a role at frequencies different from that 
given by the Bragg condition: 


oz 




(29) 


which yields 


^(*9) _ ( r(s9) 
^ki ^ki 


( 30 ) 


where, the overlap integral is obtained after approximat¬ 
ing w = d (see Eq. 0) , as 


r(s9) 

^kl 


1 


uojfi 


d 

dxdyEk ( Sx ——h x 
dx 


,d^ 

dSP 


El (31) 


Following the arguments yielding the Bragg condition 
Eq. (28), we can expand 


I) =Ed:?>e-"- 


and derive a square gradient Bragg condition 

O-TT 

Pk-Pi=m— (32) 

For backscattering into the same mode {Pk = —Pi) it 
reads 

Pk=m'^ (33) 

This looks just like the ordinary Bragg condition, but 
note, that here we have expanded the square gradient. To 
see the difference we consider a specific boundary q{z) ~ 
sin(^ 2 ;). The Fourier series of the square gradient of 
q{z) reads 


5sin(^z) 

dz 


27r 

T 


1 I 

2 + 4® 


1 ,4^ 

4® '' 


(34) 

In contrast to the Bragg condition (m = ±1), here we 
have two contributions m' = 0 and m' = 2. This means, 
that the square gradient Bragg scattering impacts the 
transmission at two disjunct frequencies. 

Plugging m = 0 into the square gradient Bragg condi¬ 
tion Eq. (33), shows that backscattering into the same 
mode (Pk = —Pi) occurs for Pk 0. In symmetric 
dielectric waveguides (as the exemplary waveguide, we 
discuss) this effect is not present because modes have a 
non-zero cut-off frequency. This means, that for small 
wave vectors there is no guided mode that could be af¬ 
fected by the m = 0 scattering. For asymmetric dielectric 
waveguides it should in principle be possible to observe 
strong backscattering due to square gradient Bragg re¬ 
flection for /3 —0. This effect has been demonstrated ex¬ 
perimentally within the statistical approach for perfectly 
conducting metallic waveguides the do not have cut-off 
frequencies m- We will show that our theoretical frame¬ 
work covers these results as well. However, contrast to 
the statistical approach [9] , which works only for bound¬ 
aries that feature peculiar statistic features our approach 
covers arbitrary shaped boundaries. In particular, we are 
not longer restricted to (pseudo) random boundaries. 

The second contribution is to = 2. It operates at half 
the wavelength of ordinary Bragg scattering. This is in 
fact a general feature. Every boundary that obeys 


q(z +A/2) = -q{z) 
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will exhibit square gradient Bragg resonances at half the 
wavelenght of the Bragg resonance. This is, because the 
square of such a q will have a periodicity A/2. 

In contrast to Bragg scattering the frequency domain 
where square gradient Bragg scattering occurs didn’t 
attract much attention, presumably because there was 
no further Bragg order expected or the square gradient 
Bragg resonance was confused with higher order Bragg 
resonances. Therefore, this reflection resonance has - to 
our knowledge - not been observed or identified experi¬ 
mentally. 


VI. COUPLED MODE EQUATIONS FOR 
ARBITRARY BOUNDARIES 

So far we have qualitatively investigated infinite peri¬ 
odic systems, that could be expanded in a Fourier series. 
In this sections it is shown how arbitrary (finite) bound¬ 
ary profiles can be treated quantitatively. We will see 
that the Bragg condition will be replaced by its continu¬ 
ous counter-part, the Fourier transform of the boundary. 
So, instead of expanding the boundary as a Fourier series, 
it will now be represented as a Fourier transformation. 
This means dropping the assumption of a periodic func¬ 
tion. 

9(5) = E ^ f (35) 

m 

To obtain the Bragg condition in the continuous case, the 
dAfc is, as in Eq. ( [26| , integrated over a small domain s: 

dAu = -i^^Y. J dfdqtW 

’■ -V-^ 

-i.O,V/3fc-/3!//3 

(36) 

The exponential function in the second integral oscillates 
and will thus be zero, as long as Pk — Pi 7 ^ P- There¬ 
fore, the result can be approximated by a normalized 
(5-Function 

dAfc « E / dPqb{P)N6{Pk - Pi-P) Jdz lifAi 

(37) 

= -^^\^d^^biPk-Pi)J^dzli^;>Ai (38) 

and in case of backscattering, where Pk = —Pi 

= -*|^E^9^(2/3fc)^dz4?)A, (39) 

The normalization N is discussed below. The coupling is 
thus proportional to qb{2Pk). This line of reasoning holds 
for both qb{P) (Bragg scattering) and qsgiP) (square gra¬ 
dient Bragg scattering). Instead of discrete conditions 


Eq. (28) and Eq. (33) we now have a continuous range 
where scattering can occur. This continuous range is 
given by the Fourier transformations qb/sg of the bound¬ 
ary and the square of the curvature of the boundary. Tak¬ 
ing the derivative of Eq. (38) yields a generalized coupled 
mode equation for arbitrary boundaries (compareEq. ([^) 


= (40) 


As before, investigating two modes, yields a set of two 
coupled equations {k = 1,1 = 2 and A: = 2, ^ = 1). Solving 
this set for contra-directional coupling {Pi = —P 2 ), yields 
two solution, one for the incoming mode Ai and one for 
the backscattered mode A 2 . As before, the reflectivity is 
given by the ratio of the incoming and the backscattered 
mode. Plugging in the solutions, they evaluate to 


p(^) _ 

^kl — 


A2(0) 

Ai(0) 


= |tanh(iVgfe(/3i-/?2)4.?^A) P (41) 


for Bragg scattering. Starting from and 7^/®^ in¬ 

stead we arrive at the result for square gradient Bragg 
scattering 


R 


(sg) 

kl 


4 ^ 2 ( 0 ) 

I/M 


\tmA,(Nqsg{Pi- P2)lifL) 

(42) 


This is the main result of this paper, and shall be dis¬ 
cussed in detail. At first, it is in perfect agreement with 
the maximum reflectivity in the periodic case, Eq. <§■ 
By comparing Eq. (§ and Eq. ( |4T| ) we can fix the nor¬ 
malization N. The result of the two approaches has to be 
the same when evaluating the coupling of two modes in 
an infinite sample over a range L or in a finite system of 
length L. Thus setting Nqb{Pb) = Pm=ii yields N = ^. 

At this point we are able to show, that our results in¬ 
clude the previous results from the statistical approach 
[9] . When assuming a perfectly electric conducting hol¬ 
low waveg uide qblj^f qsglkT^ evaluate to (see Ap¬ 
pendix D) 


L 


1 

(b),{AS) 


{dsgikl ) — (b),{SGS) 

J-'n 


1 


(43) 

(44) 


for even modes. jdP)dAS/SGS) backscattering 

length derived in [5]. In contrast to [5], we find differ¬ 
ent results for odd and even modes. However, it seems 
that the authors of [3] were unaware that they studied a 
symmetry reduced version of the system. Therefore their 
results are only valid for modes with E{z = 0) = 0, i.e., 
even modes. 

This means that we derived the very same expression 
under drastically relaxed assumptions. First, our result 
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is valid for any boundary and the reflectivity is directly 
calculated from the boundary in a straightforward way. 
A cumbersome generation procedure to generate random 
boundaries that comply with the statistical requirements 
of the statistical approach (see nans]) are not neces¬ 
sary. Second, the derivation does not depend on the 
type of boundary. Previous studies were bound to per¬ 
fectly electric conducting boundary conditions. There¬ 
fore they could not be applied to dielectric waveguides, 
where the mode crosses the boundary and penetrates the 
region outside the waveguide. 


VII. REFLECTIVITY OF A DIELECTRIC 
WAVEGUIDE 


Now the reflectivity for a dielectric waveguide can be 
calculated in a straight forward manner. By inserting 
the modes of the E-Field (see Appendix B) into Eq. (22) 


andEq. (311,4, miu 


and are obtained. The Fourier 


transformation of the boundary qt, and of the curvature 
of the boundary are calculated in the [Appendix C[ In¬ 
serting these quantities into Eq. (41) and Eq. 11^ yields 
the reflectivity. As an example, we investigate the first 
mode in the waveguide shown in Fig. In Fig. a) 
and b) the calculated reflectivity R and the transmission 
T, respectively, are displayed as a function of (vacuum) 
wave vector k. Two sharp peaks caused by Bragg scat¬ 
tering into the same mode R^i and into the next odd 
mode R^i^ can be clearly identified. The new square gra¬ 
dient Bragg reflection resonances r[\^\ and occur 

as expected at twice the value of /3. From this calcula¬ 
tions it is apparent that square gradient Bragg scattering 
is as strong as Bragg scattering. Note that these square 
gradient peaks must not be confused with conventional 
higher order Bragg peaks. In this waveguide with sinu¬ 
soidal boundaries, there are no conventional higher order 
Bragg terms, since the boundary has only a single fre¬ 
quency component. 


VIII. COMPARISON TO NUMERICAL 
RESULTS 

To test our analytical predictions numerically we use 
a commercial hnite element solver (JCMwave). The 
Maxwell equations are solved on a non-uniform 2d mesh. 
Convergence was tested and confirmed by increasing the 
finite element degree up to 7 [T3]. Fig. [^compares the 
analytical to the numerical results for the structure dis¬ 
played in Fig. We see in Fig. a) that for Bragg 

scattering into the same mode , , the calculated re¬ 

flectivity is in perfect agreement with the numerical re¬ 
sults. The situation is different for the square gradient 
Bragg scattering into the same mode The peak 

is clearly visible, but overestimated by the theory. For 
inter-mode scattering R^^^ and r[^^\ the situation is re- 


Pi (1/m) 



0 5.0x10® 1.0x10’ 1.5x10’ 2.0x10’ 

k(l/m) 

Figure 2. a) Reflectivity of first mode through the dielectric 
waveguide shown in Fig. An alytical pred iction (dashed 
black, ■ - -) calculate from Eq. ( |41| ) and Eq. ( |42[ | in comparison 
to numerical data (orange, — ). b) Analytical transmittance 
calculated from reflectivity in a) as T — 1 — R (dashed black, 

- --) compared to numerical transmission (green, ^— ). The 
difference in both are the radiation losses 1 — (i? -I- T) (black. 


versed. The reflection resonance is much stronger than 
predicted by the theory, for both scattering mechanisms. 
This is surprising, because the overlap integral in , 
Eq. (19), is obviously smaller for modes where k ^ I com¬ 
pared to k = I, where the overlap is maximal. The origin 
of the strongly enhanced inter-mode scattering has to 
be investigated more thoroughly, in a framework beyond 
two-wave coupled mode theory. 

To investigate the influence of square gradient Bragg 
scattering on the transmission, we numerically calculated 
the transmission: the result (green, ^— ) is shown in 
Fig. b). It is compared to the analytical transmis¬ 
sion gashed black, - ■ calculated from the reflectivity 
T = 1 — R. As expected, the reflection resonances appear 
as gaps in the transmission. Still, the transmission shows 
strong deviation, from i?-|-T = 1, due to radiation losses. 


To investigate the influence of radiation losses, we plot¬ 
ted the losses as 1 — (i? -I- T) (black, -) in Fig. |b). 
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Figure 3. E-Field (Re(-Ej;)) at maximum of radiation losses 
at fc = 1.4 • 10^/m (compare dip in Fig. [^b). Plane waves 
radiate out of the waveguide (dashed lines). The expected 
fc-vector from Eq. ||45|) is indicated (0 = 16.78°). 


strongly underestimated. We find that Bragg and square 
gradient Bragg scattering are of comparable strength. 

The transmission through the waveguide is even 
stronger affected by the square gradient Bragg reflections 
due to radiation losses. 

Since Bragg scattering is one of the key properties in 
corrugated waveguides. A general theory which is able 
to describe all Bragg resonances is a promising new tool 
for designing optical systems based on waveguide struc¬ 
tures. For example, in directional couplers, coupling and 
reflecting gratings are combined [5]. Both gratings have 
to have different periodicity. Using the square gradient 
Bragg scattering mechanism it should be possible to use 
the same grating for coupling and reflection purposes. 
This would considerably simplify grating structures. 
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It stems from light that is coupled out of the waveguide 
by coupling to radiation modes Pout = nk sin 0 
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where n is the refractive index of the outer material and 
Pk is, as before, the propagation constant of the guided 
mode. Eq. (451 yields (with parameters from Fig. 
an outcoupling angle of 0 = 16.78° (measured counter 
clockwise from the x-axis), which is in perfect agreement 
with the numerically calculated maximum (Fig.[^. 

The second important feature are the strong scattering 


losses at the reflection resonance. Around 25% of the 
transmission is lost due to radiation out of the waveguide. 
The effect is smaller (~ 15%) for curved Bragg scattering, 
but still visible. 


A general observation when comparing the analyti¬ 
cal to the numerical data is that the reflection into the 
same mode is overestimated, while backscattering into 
higher modes is underestimated by the generalized cou¬ 
pled mode theory. The transmittance is further reduced 
by radiation losses at the (square gradient) Bragg reso¬ 
nances. 


IX. SUMMARY 


We showed the appearance of unexpected Bragg reso¬ 
nances in dielectric waveguides with corrugated bound¬ 
aries. We generalized a statistical approach to cover in¬ 
dividual systems with arbitrary boundaries. 

The analytically calculated reflectivity is in good agree¬ 
ment with the numerical results. While the position of 
the expected resonance is predicted with high accuracy, 
the strength of the square gradient Bragg scattering is 


Appendix A: Explicit calculation of the coordinate 
transformation 


The differential operator in the new coordinates 

XV i Z] 

(x, y, z) —t (— y, z) with w{z) = d + 2aq{z) 


is derived in the following. For some function q{x, y, z) 
we obtain 


d {dx d dy d 


dy 

d_ 

dz 


q = 

q = 


1 

1 

d 

dx 

dx 

m 1 

1 

d 

dx 

dx 

1 

1 

d 

dy 

dx 

dx 

d 

dz 

dx 

dx 

d 

dz 

dx 


9 = 

w ox 


dz d 

d d 
w di 
dy d dz d 

dy dy dy dz 

dy d dz d 

dz dy dz dz 

d 


dxa dq d d 


^ 

dy 


a dq d 


d_ 

—w{z)dzdx ' dz 

—r—' V 















8 


where A, and B, label the terms for better tracking. Ap¬ 
plying the derivation twice yields 


dx"^ 

52 


dP 52 

w2 552 ' 

9y2 

/ 


5^2 


9 = 




5ai 


+ x^ 




2 ^ 

Si 2 


V 


yl 2 


a dq d d 


d 


dq d 


—w dz dx dz dz —w{z) dz dx 


AB 

\ 

52 / 

where the last term becomes 

5 . a dq d 
dz —w {z) dz dx 
.2 / o X 2 


BA 


X(J 

liP 

XU 


dq 

dz 

d'^q 


d_ 

dx 

d 


dqV d_ 

\dz) dx 
XU dq d d 


xu 

ii2 


—w \dz^J dx —wdzdzdx 

There appear different types of derivatives of q. Group¬ 
ing these terms using the reduced Laplacian defined in 
Eq. (tMl) 


V 2 = V 2 


red 


sg 


we have 


V 2 = ( 4 - 1 

^2 cr dq^ 
= —^x 


5^ 
5a; 2 
5 5 


-h 


5 5 


w dz \dx dz dz dx 


u (d^q\ ^ 5 


u- I dq\\^ d u^ (dq 

^ I ~ I + — ( — 


w \ 552 ) dx 

52 


Note that only the x coordinate was transformed, hence 
z = z and y = y. 


Appendix B: Modes of the electric field in dielectric 
waveguides 

Fig. [3] shows the first mode of the electric field in the 
waveguide. The mode numbers above are chosen to be k 
and I, to identify two specific, but not necessarily differ¬ 
ent modes. In the following the general index m is used. 
In the substrate (a; < —d(2) and cover (d/2 < x) regions, 


the m-th mode is given by the E-field perpendicular to 
the x-z-plane as (see 0) 

^ 1 -Es(fc) exp ( 7 s (i-I-D) in substrate 

— 1 i;, /;\ / d\\ 

I Ec{k) exp (— 7 c [x — j m cover 

and in the film region {—dj^ < x < d/2) odd and even 
modes are given by 


Emix,k) = Ef{k) 


cos {nfx) for odd modes 
sin(K/a;) for even modes 


with 


^2 = nV -Bl = 


= njk^ - /3^ 

^2 = n2fc2 -dl = 


■ Pv 

Here, the refractive indices refer to the substrate (ug), 
cover (rZc), and film (n/) material. The wave vector 
ricj^sk has components in direction of propagation j3m 
and in transversal direction denoted by K,c,f,s- Note, that 
there is one 5m for each mode. Consequently k and 7 de¬ 
pend on the mode number just as Em- In the case of a 
symmetric waveguide the allowed 5 m for each mode m 
can be found by numerically solving (see [S]) 


2dK / — 4 arctan 


1 - 


f — 

V"=// J 


\{^y 


= 27rm 


- 1 


where rig// = 5m/ 1 ^ 1 - The total normalization is chosen 
as 


25m 

uojy 


J dxEm = IW/m 


and yields the peak field inside the waveguide for 
odd/even modes 


^2 ^ sin(d«;/) 


5n 




- ^Iff 1 

n} - nl 7s 


- <ff 1 

n}-nl 7c 


W/r 


The amplitudes are connected via 

E){n) - niff) = Elin} - nl) = El{n) - nl) 


Appendix C: Fourier Transform qb and qsg of the 
boundary 

Let the (normalized) boundary q of the waveguide (see 
Fig. 0 be of the simplest possible form: 


q{z) = •\/ 2 sin( 5 b 2 :)boxQ (z) 
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where /?{, = ^ is the periodicity of the boundary and 
boxQ{z) is a box, or rectangular, function of length L 
constructed via the Heaviside function 0, as boxQ(z) = 
0(z) (1 — 0(z — L)). The wave vector of the boundary 
roughness /3f, should be chosen in such a way that 
are multiples of 27r to ensure a continuous function. The 
Fourier transform is normalized such that 

q{l3) = FT {q{z)) = ^ J q{z) exp (i^z) 


To calculate qb and qsg the following derivatives with 
respect to z have to be calculated 

q' = cos(,0bz)boxQ (z) + •\/2sin(/3bz) {S{z) — S{z — L)) 
q'^ = 2 ^/3bcos{,5bz)boxo (z)^ 

+ 2(^sin(/3f,z) ((5(z) - 6{z - L)) ^ 

+ 4,^6 cos(/3f,z)boxQ (z) sin(/3bz) (5(z) - 5{z - L)) 

Now the Fourier transformation of q'^ has to be calcu¬ 
lated. However, the boundary profile is chosen in such a 
way, that it ends at a; = 0, i.e., at the position of the un¬ 
perturbed boundary. Therefore sin(/3bL) = 0, and thus 
the second and the third term vanishes since 

J sin^(^f,z) {S'^{z) — (z — L)) exp{—ikz) 

= sin^(O) exp(O) — sin"^ {f3bL) exp{—ikL) 

= 0 

Using boxg (z)^ = boxg (z) and rewriting the square of 
the cosine as sum of two cosines, yields 

FT {q'^) = 2l3lFT (cos(/3bz)%ox^(z)2) 

= 2(ilFT (^(cos(O) -I- cos(2/?f,z)) boxg (z)^ 

= 2/3bFT (^boxg (z)^ -f PbFT (^cos{2f3bz)hoxQ {z)'j 

The intermediate result in Fourier representation is 

{q) = V2FT (^sin(/3bz)boxo (z)^ 

FT {q'^) = 2l3l (ft (box^(z)) + FT (cos(2/3bz)box^(z))) 

(Cl) 


With the convolution theorem 

FT{fg)=FT{f)*FT{g) 

this result can be further simplified. The (*) operator de¬ 
notes a convolution. Especially interesting for the present 
case is the convolution with a (5-function, which evaluates 
as 


f{z)*S{z-b) = f{z-b). 


Shifting a function, will add an additional phase factor to 
its Fourier transformation. Therefore instead of evaluat¬ 


ing the functions in Eq. (Cl I directly they will be shifted 
to be centered around z = 0. The Eourier transformation 
of a centered box function is given by 


FT (box^:///2(z)) = ^sinc 


(f 


and the Fourier transformation of the Cosine and Sine is 
given by 

FT (sin Pbz) = l(5(/3 - /3b) - l,5(/3 + /3b) 

FT (cos 2(3bz) = ]^5{I3- 2,5b) + ^^(;d + 2/3b) 

Applying the convolution theorem and evaluating the 5- 
functions the resulting expression reads 

-sine Q(/3-b/3b)T^^ 

FT {q'^) = ^Pb (^sinc -b ^sinc Q(^- 2/3b)T^^ 

+ ^/3b Qsinc Q(/3-b2/3b)L^^ 

Now one can neglect the last terms in both equation*, be¬ 
cause they contribute for negative frequencies only. The 
expression for q^g can be separated in two different con¬ 
tributions. The first contribution, responsible for a peak 
at ,0 = 0 is nothing but the Fourier transformation of the 
box-function. The individual shape of the boundary has 
no influence for small /3. The second term is responsi¬ 
ble for a peak at /3 = 2/3ba;, similar to the peak of qb at 
P = PbZ. 


Appendix D: hollow waveguide with perfectly 
conducting walls 

In this section we will now apply our findings to spe¬ 
cial case of hollow waveguides, with perfectly reflecting 
boundaries. We will show, that previous theoretical work 
is included in our theory. 

Assume a hollow metallic waveguide, such as a mi¬ 
crowave waveguide. It is convenient to assume a perfect 
electric conductor. That is assuming that the E-Field 
vanishes at the boundaries. The longitudinal wave vec¬ 
tor P is than a simple function of k 



P can take any number from 0 to oo. This means that 
k can have any orientation from transversal to nearly 
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longitudinal. The mode will be zero everywhere except 
inside the waveguide, between x = —d/2 and x = d/2. 
This means that the dielectric mode given in [Appendix 
[Bj is drastically simplified by 7 c = 7 s = 0 and Kf = ^. 
it is thus restricted to interior of the waveguide —d/2 < 
X < d/2: 


Now and 7^®®^ can be calculated as 

= <?^9(2/3)4r^ = (^1 - <lsg{‘2l3) for odd n 

= 9^9(2/3)4r^ = (^1 + '?«s(2/3) for even n 


Ey{x, k) 


Eoik) 


cos i^x) 
sin(^i) 


for odd modes 
for even modes 


Normalizing the power to 1 W/m yields Eq = 


After identifying |gf,(2/3)p with lT(2/3) and \psg{2j3)\'^ 
with 2S{2/3) [12], we see a surprisingly simple relation¬ 
ship between coupling coefficient and localization length 
Ln for the even modes: 


1 

rib),(AS) 
J-'n 

1 

rib),iSGS) 

J-/n 


(c<r7 
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